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ABSTRACT 

We study the evolution of the vector magnetic field and the sunspot motions observed in AR 6555 during 
1991 March 23-26. This region displays two locations of large magnetic shear that were also sites of flare 
activity. The first location produced two large (X-class) flares during the period covered by our observations. 
The second location had larger magnetic shear than the first but produced only small (M- and C-class) flares 
during our observations. We study the evolution of the photospheric magnetic field in relation to the large 
flares in the first location. These flares occurred around the same included polarity and have very similar 
characteristics (soft X-ray light curves, energies, etc.). However, the whole active region has changed substan- 
tially in the period between them. We found several characteristics of the region that appear related to the 
occurrence of these flares: (1) The flares occurred near regions of large magnetic “shear,” but not at the loca- 
tions of maximum shear or maximum field. (2) Potential field extrapolations of the observed field suggest that 
the topology changed, prior to the first of the two flares, in such a way that a null appeared in the coarse 
magnetic field, (3) This null was located close to both X-class flares and remained in that location for a few 
days while the two flares were observed. (4) The flaring region has a pattern of vector field and sunspot 
motions in which material is “squeezed” along the polarity inversion line. This pattern is very different from 
that usually associated with shearing arcades, but it is similar to that suggested previously by Fontenla and 
Davis. The vertical electric currents, inferred from the transverse field, are consistent with this pattern, (5) A 
major reconfiguration of the longitudinal field and the vertical electric currents occurred just prior to the first 
of the two flares. Both changes imply substantial variations of the magnetic structure of the region. On the 
basis of the available data we suggest that these changes made the flaring possible, and we develop a scenario 
that can explain the origin of the magnetic free- energy that was released in these flares. 

Subject headings: Sun: flares — Sun: magnetic fields — Sun: X-rays, gamma rays 


1. INTRODUCTION 

An outstanding problem in the physics of the outer layers of 
the Sun is the physical mechanism that builds up magnetic free 
energy and releases it in solar flares. The production of major 
flares is expected to be closely related to some particular 
properties of photospheric magnetic fields. Suggested proper- 
ties include particular values of the magnetic field, gradients, 
shear, sunspot motions, and magnetic field patterns that may 
indicate the probability and magnitude of an impending flare. 
The identification of the relevant properties is important for 
developing and testing theoretical models, and for forecasting 
flare activity. 

However, observations point in different directions regard- 
ing the properties of flaring active regions. For instance, it has 
been long known that, statistically, the most flare-active 
regions are the “delta” configurations where opposite polarity 
sunspots are intermingled in a common penumbra (e.g., 
Kunzel I960). Also, it is known that the most dynamic regions 
are more flare prone (e.g., see Tandberg-Hanssen & Emslie 
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1988), and more recently it has been shown that many flares 
occur in highly “ sheared ” active regions (Hagyard et al. 1984). 

All the above mentioned properties derive from statistics 
over many active regions, and sometimes a particular flaring 
active region may not show these properties. Moreover, some 
active regions having these properties produce no large flares. 
Known indicators of a flare-prone situation may not, therefore, 
fully determine the occurrence of a flare. Other active-region 
properties that are sometimes related to flares are emerging 
and canceling flux (e.g., Heyvaerts, Priest, & Rust 1977; Livi et 
al. 1989). Their roles are still controversial, and a number of 
examples and counterexamples are found in the literature. 

Each of the previously mentioned characteristics of the 
flaring active regions has influenced theories of flare buildup 
(e.g., see reviews by Svestka 1976 and Sturrock 1980): emerging 
flux is central of the model by Heyvaerts et al. (1977), in which 
a current sheet forms as a result of the action of new flux 
against preexisting flux; magnetic shear spawned a variety of 
flare scenarios based on the sheared arcade (e.g., Sakurai 1989) 
or a twisted loop (e.g., Spicer 1981). Other studies suggested 
that flares arise from the interaction between loops in adjacent 
active regions (e.g., Machado 1987). Recently two particular 
features have been related to the occurrence of major flares: 
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according to Leka et al. (1993) and de la Beaujardiere, Can- 
aeld, & Leka (1993), the locations of vertical currents derived 
from pholospheric vector magnetograms are close to flare 
kernels; and Demoulin, Henoux, & Mandrini (1992) suggest 
that magnetic “ separatrices ” are the locations where flares 
originate. 

Many of the above-mentioned characteristics and models 
appear, at first sight, as unrelated, and some are mutually 
exclusive from a theoretical standpoint. There is no compelling 
evidence that all flares, of all magnitudes and various spatial 
and temporal behavior, must share the same configuration of 
the magnetic field and identical process to build up magnetic 
free energy. Moreover, some disparate cases have been 
observed in which it seems that no single model may apply to 
all stages of the observed flares. Therefore, we must keep an 
open mind and accept that similar basic plasma processes for 
energy storage and release may occur that lead to various 
kinds of flares in several different magnetic configurations 
(Gaizauskas 1989). We need to examine many cases in detail 
for establishing the ways an active region may develop to a 
flare-prone situation, and to extract the relevant properties of 
the magnetic field. 

In this paper we will not try to single out a “cause” or a 
particular scenario for all flares. Rather, our goals in this paper 
are to analyze (1) how the interpretation of a set of observa- 
tions may indicate magnetic free-energy buildup, and (2) how 
the observations may indicate the likelihood of occurrence of a 
particular type of flare. Also, we compare the magnetic struc- 
ture of flaring and nonflaring parts of an active region to assess 
the role of various properties related to flaring. In particular 
we study the shearing of the field, and the large vertical electric 
currents, that did not lead to flaring. 

We study observations of the vector magnetic field of the 
active region AR 6555 for the period 1991 March 23-26, per- 
formed at MSFC. These data are supplemented by white-light 
and sunspot structure and position from the Debrecen Obser- 
vatory. Earlier study of this region, during the same time 
period, by Ambastha, Hagyard, & West (1993) concentrated on 
the characteristics of the observed “ shear” through the overall 
region. Our present study uses similar data but concentrates 
on different aspects concerning the field topology, the develop- 
ment of electric currents, and the development of magnetic free 
energy that may power the flares in the active region. 

In the next sections we will describe the basic structure and 
parameters of the active region and flaring we study. Then we 
will describe the observed sunspot motions, and the evolution 
of the coarse structure of the potential field that is inferred 
from the vector magnetograms. In subsequent sections we 
study the observed departures of the field from the potential 
case. From these departures, and the sunspot motions, we infer 
how, in the flaring region, magnetic free energy appears to 
build up in a very specific way. Finally we discuss our interpre- 
tation of the observations and their implications for flare 
models and for flare forecasting schemes. 

2. THE ACTIVE REGION 

AR 6555 passed across the visible hemisphere from 1991 
March 17 to 31 at a solar latitude around 25°-28° South. Our 
study is centered around the central meridian transit from E20 
to W20, and thus is not seriously affected by projection effects. 
The white-light appearance of the region is shown in Figure 1 
(Plates 12-13), which shows the main sunspots of the region 
labeled for future reference. 


2.1. The Flaring 

In the days preceding our observations flares of M- and 
C-class occurred, and also an X-9.1 event occurred on March 
22 at 22.45 UT. Most of these flares were located in the leading 
complex around P3 and N6, and some of them were quite 
extended (see Ambastha et al. 1993). 

During the period of our observations, in contrast, no large 
flares (X-class) occurred near P3 and N6, even when large 
shear, vertical currents and field gradients existed at these loca- 
tions as we show below. The only large flares reported were 
two X-class flares located near PI: the first on March 25 at 
08:10 UT (X5.3), and the second on March 26 at 20:35 UT 
(X4.7), In this paper we study the evolution of the magnetic 
field in relation to the occurrence of these two large flares and 
the absence of large flares near P3 and N6. 

2.2. The Observed Magnetic Field 

The evolution of the longitudinal field structure of the region 
during the observed period is depicted by Figure 2a (Plate 14). 
The largest magnitude of the negative longitudinal field is 
reported at about 2 100—2200 G by Mount Wilson (the MSFC 
data is suitably corrected to match these data, as described by 
Ambastha et al. 1993), and corresponds to the large umbra, Nl, 
of the main trailing sunspot group. The largest positive field, 
reaches only about 1800 G at one of the small leading spots, 
P5. These spots are well separated from the main trailing 
group by an area of intermingled polarities. In Table 1 we list, 
for reference, the maximum values of the field for each of the 
spots identified in Figure 1. A large flux imbalance is evident: 
the accumulated positive (leading polarity) flux can only 
balance 15%-18% of the accumulated negative (trailing 
polarity) flux. This result cannot be accounted for by any rea- 
sonable calibration offset. Consequently, we expect that most 
of the negative flux closes outside the field of view of our 
magnetograms. There is a neighboring active region farther to 
the east, another small active region to the north, and a few 
small positive spots to the west, but they are not enough to 
compensate for the missing return positive flux. Therefore, we 
conclude that much of the return flux should consist of small 
magnitude longitudinal magnetic field that is extensively 
spread over a large area, probably through the enhanced 
network that extends on the east in the Ca ii images. 

In Figure 2b we present a saturated view of the longitudinal 
magnetic field (in which all fields larger than +200 G are 
displayed white and all fields smaller than — 200 G are black), 
in order to show the weaker fields and the main polarity inver- 
sion line (hereafter MPIL). This line is located to the east of the 

TABLE 1 

The Sunspots’ Maximum Field on March 23 

Longitudinal Field Denomination in Ambastha 


Spot (G) et al. 1993 (Fig. 2) 


PI 1126 PI 

Nl -2155 FI 

N2 -1965 F2 

N4 -1266 F4 

P2 518 P2 

N3 -1298 

P3 1455 P3 

N5 -1537 F5 

P4 813 P4 

N6 -1454 F7 

P5 1832 P5 
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main leading spots, P5 and P4, at the region of relatively low 
fields (the other main leading spot, P6, is often at the edge of 
the field of view). The figure also show the complex positive 
polarity, P3, that invades the negatively dominated area 
(formed by N3, N4, N5, and N6), the steep field gradients with 
surrounding negative field small spots, and the included spots, 
PI and P2, that surround the main leading spot group, N1 and 
N2. As shown in the previous figures, these included spots 
share the same penumbra of the main negative spots. The most 
significant included polarity related to the flares observed in 
this period, IPl, is associated with the spot PI. The included 
polarity IPl has moderate field of only -- 1 100 G (at the begin- 
ning of the period). The large-scale pattern of the field is domi- 
nantly negative polarity of the west of the MPIL, and 
essentially positive to the east. This is the expected pattern in 
this cycle and hemisphere, but it is somewhat atypical in that 
the trailing flux of AR 6555 is highly concentrated while the 
leading flux is very spread out. A diagram showing how the 
coarse large-scale configuration may be described is shown in 
Figure 3. 

Figure 2 shows the basic evolution of the structure of the 
longitudinal field during these days. The strip of positive 
polarity, IPl, as seen on March 23 and 24, broke in two parts 
(IPla and IPlb) between the last observation on March 24 and 
the first on March 25 (at 13:57 UT). The break was at the 
southeastern-most point, and this location changed from being 
part of IPl (positive polarity) into being part of the main trail- 
ing sunspot (neptive polarity). This change cannot be 
explained by projection effects, because the region is not far 
from centra] meridian (and close to disk center), and because of 
the direction of the transverse field (almost perpendicular to 
the region displacement on the disk between the two days). The 
basic structure in this part of the active region remained 
unchanged through March 25 and 26, but the spot PI associ- 
ated with one of these fragments (IPla) drifted away relative to 
the main negative spots until it finally detached from the 
common penumbrae. 

The first large flare near this location (on March 25) 
occurred close to the time of the splitting of IPl; a precise 
timing is not possible because these events occurred during the 



Fig. 3.— Scheme showing a way in which the overall large-scale field can be 
represented by a tilted submerged dipole. The vertical negative flux is very 
concentrated and leads to a large field in a relatively small region, while the 
positive flux is very spread and corresponds to weak fields. Note that the field 
to the east is mostly horizontal. 


night at the MSFC the first observation on March 25 is about 
5:47 hr after the flare). The longitudinal field undergoes slow 
fading after this rupture; however, the transverse field changed 
more significantly. This change is such that the magnetic shear 
at IPla decreased (as did the electric current; see below). 
However, the second large flare occurred on March 26 at basi- 
cally the same location and with very similar energy and tem- 
poral behavior. This flare and the vector magnetic fields were 
observed at MSFC with good temporal resolution (see 
Ambastha et al. 1993). 

3. SUNSPOT MOTIONS 

We derived the proper motions of umbrae in AR 6555 
during 1991 March 21-26 from 203 photographic observa- 
tions, taken at the Debrecen Heliophysical Observatory and its 
Gyula Observing Station. Details of the method for reducing 
the measurements are given by Kalman (1980), and Bumba et 
al. (1993). The heliographic coordinates of the various umbrae, 
computed in the Carrington System, show a strong backward 
drift caused by differential motion at the latitude of the sunspot 
group. Proper motions of the entire group are difficult to assess 
because of uncertainties in the differential rotation and in dif- 
ferences between the rotation rates of different solar features. 
For this reason we consider here only the motions relative to 
the largest field and most stable umbra Nl. Since the individ- 
ual measurements have a scatter of about 0.15 heliographic 
degrees, mainly due to atmospheric seeing effects, and there are 
more than 200 measurements of each long-lived umbra in the 
course of 6 days, some averaging and smoothing was per- 
formed to show the sunspot trajectories. We fitted cubic splines 
to the mean positions of the umbrae at 12:00 UT on each 
observing day. 

Figure 4a gives the trajectories of the individual umbrae in 
AR 6555 that were clearly observed for at least 2 days. In this 
paper we designate the umbrae in a way resembling the one 
used in Ambastha et al. (1993), namely, the leading polarities 
(positive) are indicated by a prefix P and the following 
(negative) by a prefix N. Arrowheads represent the position of 
a given umbra, and its direction of motion, at 12:00 UT every 
day. The contours in Figure 4b correspond to the umbrae and 
penumbrae, with heliographic coordinates corrected for differ- 
ential rotation. These contours show the morphological 
changes in the active region during 1991 March 21-26. 

AT 6555 was dominated by the large complex Nl of irregu- 
larly shaped, old following-polarity umbra dissected by several 
light bridges- This complex moved eastward slowly with 
respect to the local standard rotation as is characteristic of old 
spots. (This eastward motion is not shown in the figure, where 
all motions are referred to those of the Nl complex.) The most 
significant changes occurred in the eastern part of the group, 
where a new pair of polarities (P1-N2) appeared on March 21 
in an almost N-S alignment. This new dipole moves to the west 
as is usual for new spots, and its northernmost part, N2, moves 
toward the similar polarity old umbra Nl. The southernmost 
part of the new dipole, the included polarity PI, emerged 
within the southeast penumbra of the old spot Nl. This spot 
grew with a rapid westward motion in the first days, and later 
moved southward, until it finally detached completely from Nl 
(on March 26). The other included spot, P2 located to the 
northwest of N2, also displayed westward motion but of 
smaller magnitude. 

On the western side of Nl spot, N5 and P3 moved signifi- 
cantly: umbra N5 with rapid westward direction; P3 with 
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ISO IS5 190 195 
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Fig. 4b 

Pig 4 — Scheme showing the sunspot motions for the period between 1991 March 21-26. The arrowheads show the positions of the umbrae at 12:00 UT on each 
day, and the directions of the arrows indicate the direction of motion. The motions are shown relative to the position of the largest field umbra, Nl, that is indicated 
with a cross. Panel (a) shows the umbrae’s positions, and panel B shows the intensity contour shapes reduced to the heliographic Carrington System. Many of the 
small umbrae whose displacement is shown here were not labeled before, and one of these umbrae to the South of N 1 displaces crossing the path of the spot PI . 


slower motions. The relative motions of these spots tended to 
bring them together. During March 23-24, the penumbra of 
N5 joined that of P3; N5 slowed down and turned southward 
until both spots moved practically parallel to each other. 

The general picture of observed sunspot motions resembles 
the pattern of a hydrodynamic flow around an obstacle. The 
role of an obstacle is played by the old, disintegrating large 
umbra Nl. The emerging material that carries the new mag- 
netic flux, shown by N2 and PI, flows around the “obstacle” 
in its westward motion. Other smaller negative polarity 
umbrae to the south of Nl also take part in this westward flow. 
As Figure 4a shows, the trajectory of the spot N7 crosses that 
of PI, but N 7 moves to the west of PI, and the two spots never 
approach each other. 


The full detail of the sunspot motions is given here because it 
is a critical test for understanding the magnetic field evolution 
and the scenario that we address in the discussion. 

4. THE INFERRED COARSE POTENTIAL FIELD 

We computed the coarse potential field that corresponds to 
the observed vector fields by averaging the measurements over 
4x4 pixels (i.e., about 11" x 11"). Averages made over 2x2 
pixels show no significant differences for our present dis- 
cussion. Since our fields of view differ for all the dates con- 
sidered, we extract a sub-field-of-view that was always 
observed and contains all the relevant features. For this sub- 
field we carried the averaging of the measured values of all 
components of the vector field and assigned the resulting value 
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to the center of the 4x4 pixels. This procedure is somewhat 
crude but it is adequate for our purpose of finding the coarse 
structure of the potential field. 

Details of our procedure for calculating the potential field 
are described in Appendix A. We describe the potential field by 
assigning a number of source dipoles to a plane at 8 pixels 
depth beneath the photosphere (twice the separation between 
the used measurements of the field). The number of dipoles 
equals the number of points to be matched. In contrast to 
other methods that consider only vertical dipoles (e.g., Demou- 
lin et al. 1992), we consider dipoles with three arbitrary com- 
ponents of the magnetic moment (i.e., dipoles with arbitrary 
strength and tilt). We find the strength of all three components 
of the dipoles in such a way that their field matches the 
observed longitudinal field exactly and matches the transverse 
field approximately. (The approximation of the transverse field 
is in the least square sense.) The horizontal components of the 
field contain important information on how the field lines 
close. They supplement the lack of constraints on the normal 
components of the field at the sides and the top of the domain 
where we compute the potential field. (Keep in mind that the 
potential field is uniquely defined only when the normal com- 
ponents of the field are given for a closed domain where no 
currents exist.) If the observed horizontal components are not 
used, and one resorts to purely vertical dipoles, the result for 
our region in which a large flux imbalance occurs would be 
unrealistic. The global structure would depart too much from 
the leading positive and trailing negative polarities that are 
apparent in large field-of-view magnetograms. Our computed 
distribution of dipoles contains only a weak vertical com- 
ponent and a much larger horizontal component of the dipole 
moment. 

The resulting potential field lines are complicated because of 
the mixing of polarities and the complicated distribution of the 
photospheric field. Since a limited field of view is used, and 
large nonpotential behavior is observed at some locations, the 
potential field calculated here represents only the one that is 
closest to the observations. The real coronal magnetic field is 
expected to differ from the calculated one because of edge 
effects, and because of electric currents that are knawn to flow 
across the photosphere and close within the upper layers. 

However, we found an interesting property of the potential 
field (Table 2). Initially the coarse structure of the potential 
field evolves only in magnitude but not in topology. No null 
points of the potential field occur above the photosphere, in 
the neighborhood of I PI, for the data corresponding to March 
23. But, on March 24 (several hours before the flare) a local 
minimum of the magnitude of the field or quasi-null, B ~ 7 G, 
appears above the location in the photosphere were the 
polarity IPl later broke up. Then, on the 25th (^6 hr after the 
first flare) the local minimum evolved into a full null of the 
field, having zero absolute value, and shifted to the southwest. 
For the potential field of the 26th, the null remained more or 

TABLE 2 

The Nulls in the Region above IPl 


Date Time Height 

(1991) (UT) (km) 


March 23 18:03 

March 24 16:41 2500 7 G 

March 25 13:57 2000 0 

March 26 13.59 1800 0 


less at the same location but shifted slightly further to the west. 
These data show that a large qualitative change in the coarse 
potential field topology took place between March 24 and 25, 
in coincidence with the onset of the flare activity at the location 
of IP la. 

The qualitative change is also apparent in some of the field 
lines (not shown) that can be computed from the potential field. 
The changes of the potential-field structure after March 25 are 
not very significant. Instead, the changes between March 24 
and 25 are substantial and correspond to a split of the eastern 
field lines from the eastward directions they have on March 24 
into northward and southward lines. These changes, and the 
appearance of the potential field null, are consequences of the 
changes that the photospheric field was undergoing, particu- 
larly when IPl split into IPla and IPlb. The change in the 
structure of the potential field seems spatially and temporally 
related to the beginning of major flare activity at IPla, 
although the theoretical reasons are not clear (but see the dis- 
cussion in § 7). Furthermore, the fundamental changes that the 
active-region structure was undergoing are hinted already by 
the appearance of a quasi-null late on March 24 (many hours 
before flaring starts). These changes were not apparent in the 
raw magnetograms on March 24, but were clear in the data 
corresponding to the following day. 

5. DEPARTURES FROM THE POTENTIAL FIELD 

The departures of the observed field from a potential one 
can be characterized in two ways: by the angles between the 
observed and potential transverse field, or by the inferred verti- 
cal electric currents. The angle is relatively unaffected by cali- 
bration of the measurements in the regions of weak 
longitudinal field. The vertical currents are, however, strongly 
affected by calibration of the transverse field. Both the angle 
and the currents are dependent on the arbitrary criteria used 
for resolving the 180^^ ambiguity of the transverse field 
azimuth. 

In addition, any comparisons rely heavily on the somewhat 
arbitrarily selected potential field (see Appendix 1), Our 
analysis uses the potential field that was calculated as 
described in the previous section. However, we also compared 
our results for the departures with those corresponding to the 
potential field obtained following a previous scheme (Teuber, 
Tandberg-Hanssen, & Hagyard 1977). We found that the qual- 
itative behavior of the shear angle does not differ much (except 
at a few locations). More significant differences occur in the 
magnitude of the transverse field; our method gives generally 
smaller values of the transverse field that are closer to the 
observed. Since nonpotential behavior is clear in either case, it 
is hard to make a strong case for the appropriateness of any 
scheme for computing the potential field. The potential field 
derived in this paper gives the photospheric vector field that is 
closest to the observed field. Therefore, the departures of the 
observed from the potential field are minimized in our scheme. 

The angle between the observed field and the calculated 
potential field, near the inversion lines of the longitudinal com- 
ponent, has been characterized as the “ shear ” of the observed 
magnetic field. This idea was extended to locations removed 
from the inversion line by Ambastha et al. (1993) for this par- 
ticular active region. 

In Figure 5 we show contours of the vertical electric current 
inferred from the MSFC vector field measurements for the four 
consecutive days. The vertical electric current density, is 
determined as indicated in Appendix B. Comparing the verti- 
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Fig. 5.— Contour levels for the electric current density corresponding to the observations shown in Fig. 2. The solid lines correspond to upflowing currents and 
the short-dash lines correspond to downflowing current. The long-dash lines correspond to the polarity inversion lines. This map was obtained by smoothing the 
original data by averaging over 2x2 pixel ( ^ 4000 x 4000 km^) areas. The current contour levels start at ± 4 mA m " ^ and are separated by 8 m A m 


cal current with the “ magnetic shear ” shown by Ambastha et 
al. (1993) in their Figure 3 it is dear that contours of maximum 
currents occur at some of the edges of contours of maximum 
shear. 

The temporal changes observed in the current’s pattern are 
shown in Figure 5. These changes are much more drastic than 
those displayed by the shear, but the directions of their changes 
coincide. These indicate that although currents are localized, 
they affect large areas. On March 23 there are several kernels 
of large currents near the polarity inversion line between IPl 
and Nl, hereafter current system A. The most extended kernels 
correspond to positive values (upflowing currents); only part of 
the negative return currents is visible. There is another current 
system at the line that half encloses the positive polarity 
around P3, hereafter current system B, and there is a third one 
on both sides of the polarity inversion line between P2 and N2, 
hereafter current system C. Systems B and C are much more 
balanced than A, and they display more or less as much 
upflowing as downflowing current. The quantitative behavior 
of the currents is shown in Table 3 for the three main current 
systems. 

On March 24 system A fractured in several parts, decreased 
from the values discussed above, and almost half-encircled the 
spot N 1 . System B gained strength and became more compact 
as its sunspots approached each other. System C decreased 
somewhat in intensity and became more fragmented. At this 


time the strongest current system in the entire active region is 
B, while system A is far less intense. 

On March 25, after the first of the two X-class flares, system 
A changed drastically. This change coincides with the 
restructing of the longitudinal field and with the motions of the 
spot PI across the penumbra of Nl. Since all these changes 
were observed simultaneously one cannot draw any conclusion 
regarding a cause and effect relationship. The changes were 
most likely the result of photospheric motions and deeper 
seated forces acting on the photospheric plasma and magnetic 
field. The structural change of system A replaces the previous 
pattern with a pair of oppositely directed semicircles of strong 
currents that are more or less in balance and another pair of 
oppositely directed currents displaced to the west. The changes 
in the other current systems are much less impressive. The 


TABLE 3 

The Main Electric Currents on 1991 March 23 


System 

Total 

Upflowing 

Current 

(A) 

Total 

Downflowing 

Current 

(A) 

Maximum 
Upflowing 
Current Density 
(mA m~2) 

Maximum 
Downflowing 
Current Density 
(mA m“2) 

A 

1.2 X 10*^ 

3.3 X 10’^ 

40 

33 

B 

2.1 X 10*2 

1.7 X 10'* 

22 

33 

C 

2.4 X 10*2 

2.5 X 10'* 

17 

21 
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northern components of system B diminished and the southern 
components intensified; system C practically dissolved, 
showing only a pair of small kernels. At this time system A was 
as strong as system B. 

On March 26, before the second X-class flare, current system 
A had simplified into two strip currents : one upflowing in the 
north, the other downflowing in the south. These strips were 
located between N1 and the now fully detached PI. The inten- 
sity of this current system had severely decreased, consistently 
with the decrease in the “shear index ” at this location. Current 
system B continued its trend of evolution; the northernmost 
portions decreased and fragmented, while the southernmost 
portion enhanced and organized- At this time the remnants of 
system C were barely visible above the noise. 

Comparing the sunspot motions and the vertical currents, it 
is evident that they coincide, and are well related to the 
observed transverse field pattern (e.g., see Fig, 5 in Ambastha et 
al, 1993). All these quantities indicate that the field in the 
region between I PI, and later IP la, and N1 was stressed by 
horizontal flows in the southwest and west directions. The 
flows have maximum velocity near the polarity inversion line 
and also produced some rotation of magnetic features. The 
motion pattern was quite different from that usually assumed 
in numerical simulations of shearing arcades. In Figure 6 we 
compare (looking down toward the photosphere) the usual 
flow and sheared field pattern with that proposed by Fontenla 
& Davis (1991), The conventional shear pattern has zero 
velocity at the polarity inversion line, a maximum absolute 
value on each side of this line, and a change of sign across the 
polarity inversion line. This pattern produces field distortions 
and electric currents that are very different from those 
observed near IPl and IP la. The pattern of vertical current 
corresponding to the usual shear flow also changes sign at the 
longitudinal field polarity-inversion-line and also reaches a 
maximum absolute value on either side. 

The flow, transverse field, and vertical current patterns that 
Fontenla & Davis (1991) propose are quite different: a 
maximum of the current density at the inversion line, at the 
location of maximum velocity, and much weaker and spread- 
out return currents at both sides. This pattern is much closer to 
the observed, although the maximum current and horizontal 
velocity are not exactly at the longitudinal field inversion line, 
but are offset toward IPl (later IP la). The return currents are 
not observed, probably because they are masked by noise due 
to their more spread-out distribution. 


polarity inversion 



6. MAGNETIC FREE-ENERGY BUILDUP 

The buildup of free-magnetic energy in the upper atmo- 
spheric layers can be studied by computing the work done by 
the Lorentz force at the photospheric layers. This quantity 
corresponds to the Poynting vector at the photosphere (see 
Krall & Trivelpiece 1973) and gives an upper limit to the rate 
of magnetic free-energy buildup. The part that is available for 
flaring is stored above the photosphere in the form of field- 
aligned currents, without being dissipated or dispersed. 

In the following, for simplicity, we will neglect the angle of 
the line of sight respect to the vertical. 

The power upflow can be computed from the magnetic field 
and the electric current that give rise to the Lorentz force, and 
the plasma velocity. From our observations we can only 
deduce horizontal velocities, and vertical electric currents, 
Jj. Vertical velocities and horizontal currents are undeter- 
mined by our observations. Consequently the full power sup- 
plied cannot be found. From our data, however, we can obtain 
one of the components of this power. 

The vertical component of the current, together with the 
horizontal component of the field, give part of the horizontal 
component of the Lorentz force, — Another contri- 

bution to the horizontal force may come from horizontal cur- 
rents and vertical fields (except at the polarity inversion line 
where the vertical component of the field vanishes). Thus, the 
horizontal Lorentz force plotted in Figure 7 may not represent 
the total horizontal Lorentz force at locations remote from the 
polarity inversion line (if horizontal currents are present). 
Figure 7 shows significant horizontal forces very close to the 
polarity inversion lines; thus they practically correspond to the 
total horizontal lorentz forces there. Note their different align- 
ment with respect to the polarity inversion line in current 
systems A and B. On March 23 and 24 the forces at system C 
are smaller and have a similar pattern to those at B; then they 
practically vanish on March 25 as the remainder shifts away 
from the inversion line on March 26. The differences between 
the forces at systems A and B are closely related to the pre- 
viously discussed differences in the current patterns. 

The forces at system B, on the northern side on March 23 
and 24, and on the southern side on March 26, point away 
from the polarity inversion line. Their largest magnitude is 
about 10 ^ N m there are opposing forces of comparable 
magnitude at each side of the inversion line. Since the observed 
horizontal motions press the material toward the inversion 


polarity inversion line 



V 


(a) (b) 

Fig. 6.— Scheme comparing (a) the usually adopted shear of simple arcades across a longitudinal field inversion line and (fr) the form of shear that was proposed 
by Fontenla & Davis (1991). This figure shows a view down toward the photosphere and illustrates the distortion that a strai^t field line would experience under the 
shown velocity patterns. We also indicate the electric currents that correspond to the distortions of the field. 
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Fig. 7.— Arrows showing the Lorentz force component that results from the longitudinal (practically vertical) electric currents in Fig. 6, and the corresponding 
transverse (practically horizontal) magnetic field. The contours correspond to the longitudinal field levels of ± ICIO and ± 1000 G. The forces computed correspond 
precisely to the total transverse component of the Lorentz force only at the locations where the longitudinal field vanishes. Wherever a longitudinal field component 
is present there may be another component to the transverse Lorentz force due to the product of the longitudinal field with the transverse electric current (that is 
unknown). The longitudinal component of the Lorentz force is unknown^ and is due to the product of the transverse components of the field and the current. 


line, we find that work is done by mechanical over electromag- 
netic forces, and therefore there is feeding of mechanical into 
electromagnetic energy. This corresponds to the observed 
buildup of vertical currents of opposite sign at each side of the 
inversion line, much like the usual shearing of arcades shown 
in Figure 6a. Such forces in system B are observed to increase 
in later data. They are rather strong on the northern side of B 
on March 24 and on the southern side on March 26. However, 
no large flares are observed near this current system after 
March 23 ; only C-class flares were reported. 

The behavior of system A is very different, as a consequence 
of its different current pattern. The forces here have an impor- 
tant component directed along the polarity inversion line. On 
March 23-24 most of the forces are directed northward, reach- 
ing a value of about 1.9 x 10~^Nm~^ (on the 23), while a few 
pixels have a southward direction. These forces have a signifi- 
cant component acting against the observed southwesterly 
motions. Despite drastic changes in the structure of the longi- 
tudinal field between March 24 and 25 some of these character- 
istics of the force remain. In current system A there is transfer 
of mechanical into electromagnetic energy, but this transfer is 
very different from that in the usual shear pattern shown in 
Figure 6a. In the period March 23-24 the situation in current 
system A resembles that depicted in Figure 6b. Then, on March 
25, the forces are mostly transverse to the southernmost part of 
the inversion line and diverge. While the earlier behavior per- 


sists in the northern part of system A, the southern part evolves 
toward the same behavior as system B. 

We find a maximum value for the magnitude of the forces in 
Figure 7 of about 1/10 of the gravitational force at the layers 
where our magnetic field is measured (height ~ 250 km corre- 
sponding to p 6 X 10"^ kg m”^). This horizontal force is 
large compared with inertial forces, but is small compared with 
gravity. The energy transfer corresponds, therefore, to horizon- 
tal pressure gradients that are smaller than the gravitational 
pressure gradient. Assuming a scale height /i 10^ m for the 
forces (similar to the density scale height), we find from the 
observed component of the Lorentz forces in system A that the 
total horizontal push on the plasma is about 2 x 10^^ N. 

The power “pumped” into magnetic free-energy of the 
plasma contained in a layer of depth h at current system A 
depends on the velocity component parallel to the Lorentz 
force (but in the opposite direction). Assuming an upper bound 
for this velocity of 10^ m s" \ this power is about 2 x 10^^ W, 
that is, 1.7 X 10^** J day"^ The energy buildup over a day is 
marginally enough to supply the combined energy released in 
both flares (viz., ~10^‘ ergs). Furthermore, purely horizontal 
plasma motions of such magnitude would result in displace- 
ments of about 8.7 x 10^ m day~V The observed umbra’s 
motions are far less, about 7 x 10^ m day " ^ (corresponding to 
80 m s"^). It is still possible that large plasma motions occur 
outside the umbrae, but it is unlikely that they could exceed 
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10^ m s~K Thus the observed magnetic free-energy buildup 
(based on Lorentz forces due to the observed vertical electric 
currents) is insufficient to power the observed flares. The dis- 
placements along the Lorentz forces in current system B are 
much smaller, and therefore the observed energy buildup is 
much smaller than that in system A. 

The rate of observed energy buildup was largest around 
March 23 at 18:03 UT, decreased somewhat by March 24 at 
16:41 UT (but remained high), and become very small on 
March 26 at 13:59 UT. This buildup appeared related to the 
source that powered the first flare on March 25 at 08:10; there 
is no comparable buildup for the second flare on March 26 at 
20:35. Thus, the observed buildup that precedes flaring in 
current system A cannot be used to forecast later activity. The 
second flare may result from energy stored in the layers at 
and/or above the photosphere until it was released 1 or 2 days 
later. The energy for the first flare can also be inferred from the 
observed large vertical electric currents, but the currents 
decayed so much by the time of the second flare that they 
cannot be considered consistent indicators of flare probability. 

It is surprising that the high energy buildup in the active 
region, although insufficient for a major flare, occurs precisely 
at the location where the flares erupt, and occurs Just before 
the first flare. The observed buildup may be Just the residual of 
a much more important energy storage occurring at deeper 
layers, or through horizontal currents that cannot be observed 
with present instrumentation. 

The observations indicate that material emerges from deeper 
layers to the east of the spot N 1 ; it flows around the south side 
of spot N1 and shears the field. Then the displacements stop; 
consequently the material must disappear by submerging 
somewhere to the southwest of spot Nl. The previous figure 
for energy buildup corresponds to stresses on the horizontal 
component (or more accurately the transverse component) of 
the field that produce a shearing of the same type as the 
observed transverse magnetic field pattern near the polarity 
inversion line. However, the emergence and submergence of 
material, if forced by pressure and gravity, respectively, distort 
the field in a way that corresponds to a buildup of horizontal 
currents. The work done in such distortions may be substan- 
tially larger than the one that we have determined above. 

7. DISCUSSION 

The magnetic field structures of whole flaring active regions 
often display large imbalance of the longitudinal field. In the 
present case this structure can be broadly described by an 
inclined dipole (Fig. 3) and does not change much in the 4 days 
we study. The significant changes that can be flare related 
occur in included polarities with a spatial scale of about T or 
smaller. Large localized electric currents are observed near the 
included polarities. The observed total electric current is 
unbalanced in the earlier observations. 

In the first part of this section we discuss which character- 
istics of these changes can be used to forecast flares. In the 
second part we discuss some scenarios for the magnetic free- 
energy buildup that may account for flares. Finally, we study 
the amount of energy buildup indicated by the observations, 
and its implications for the buildup scenarios. 

7.1. Flare Forecasts 

The similarity between the X-ray light curve for the first 
(March 25 at 08:10 UT) and second large flare (March 26 at 
20:35 UT) suggests very similar energy releases for the two 


flares. However, the magnetic configuration, magnetic shear, 
and electric currents are conspicuously different in both cases. 
This can only be explained by assuming that the observed 
photospheric vertical currents are not those directly 
responsible for the flaring. 

Our data show that flare prediction based only on 
“magnetic shear” does not predict the time and location of 
large flares. The first of two large flares occurred not at the 
highest shear location but at a secondary location of decaying 
shear. The second flare, of similar yield, occurs in the same 
region at a time of even lower shear. These considerations are 
also valid for the magnitude of the vertical electric currents. 
The largest currents were at locations without significant 
flaring, and the current was low when the second flare 
occurred. 

We find a remarkable change in the longitudinal field and 
the current pattern near the time and location of the first 
X-class flare. At this time, the longitudinal field changes its 
topology and a null of the potential field appears. The changes 
are related to motions of secondary sunspots. Also at this time, 
the electric current changes by decreasing and shifting away 
from the polarity inversion line and by becoming more bal- 
anced (see Table 3). Changes in the field and vertical currents 
are absent in the second flare, but instead sudden sunspot 
motions are observed. 

We suggest that the existence of large vertical currents, 
emerging flux and sunspot motions, are all interrelated indica- 
tors of energy transfer from subphotospheric layers. However, 
the present observations are insufficient for assessing the 
amount of energy being transferred. To gain further insight we 
need to complement the observations with indicators of a ver- 
tical plasma velocity and horizontal electric current. 

The existence of energy transfer does not determine that a 
large flare will occur, because the energy can Just be dissipated 
or disperse. We suggest that another condition is required for 
the magnetic-energy storage, at least for some types of flares : 
unbalanced vertical current must exist near the polarity inver- 
sion line. 

Our observations suggest that a flare is triggered when a 
magnetic null appears above the photosphere. In our flares, the 
apparition of the null of the coarse potential field is accompa- 
nied by sudden a change in the vertical currents that implies 
related changes in the closing currents above the photosphere. 
These changes only indicate the onset of flaring activity, and 
the full release of the stored energy may occur throughout 
several flare events. 

The existence of the magnetic null near the flaring location is 
very important. As shown by Fontenla (1993), in regions near a 
null the plasma-beta is near or above unity so that large elec- 
tric currents can be easily raised, such as the ones observed. 
These motions close to the null would lead to localized sheets 
of very large current density in which rapid dissipation may 
occur. 

7.2. Energy Buildup Scenarios 

Leka et al. (1993), and de la Beaujardiere et al. (1993) con- 
sidered flaring in a magnetic configuration similar to that of 
the flare models related to prominence eruption (e.g., Kuin & 
Martens 1986). This is clearly not the case in our flares because 
(1) the coarse structure of the field is quite different from a 
sheared arcade, (2) the vertical electric current is unbalanced 
and close to the polarity inversion line instead of a pair of 
opposite currents on both sides of the line, and (3) there is no 
filament signature in the Ha and Ca ii data we examined. 
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The part of the active region with the strongest similarity to a 
sheared arcade was current system B (see Fig. 5) which had no 
large flares. The most conspicuous surging in Ha (see 
Ambastha et al. 1993) was current system C, which had no 
large flares. 

Kurokawa (1987) reported Ha observations of sheared fields 
in two active regions: one with very large flares, similar to our 
system A; and another without large flares, similar to our 
system B. Kurokawa suggests that there are two types of 
sheared active regions, one associated with the “emergence of 
twisted magnetic flux ropes,” and another resulting from the 
‘'collision between two sunspots of opposite polarities,” as in 
our systems A and B, respectively. The emergence of material 
and the sheared field pattern in our region resemble those 
shown by Kurokawa. However, we do not describe the situ- 
ation in this way because for the observations we can only 
ascertain that the field is sheared but not that there are inter- 
nally twisted flux ropes. 

In our flares, at current system A, the magnetic configu- 
ration has a structure in which the field is distorted sideways as 
proposed by Fontenla & Davis (1991) for flare production. The 
observed horizontal motions, magnetic shear, and electric cur- 
rents, and the derived Lorentz forces all coincide with the sce- 
nario of Fontenla & Davis (see sketch in Fig. 6b). The flow 
pattern, on March 25-26, depicted in Figure 8 (see also § 3), 
results from the emergence of plasma and negative field east of 
the spot Nl, and from the migration of the previous positive 
field toward the southwest. These motions are accompanied by 
shearing of the transverse magnetic field and by changes in the 
longitudinal field, 

7.3. The Rate of Energy Buildup 

The observations reveal important Lorentz forces at the 
layers where the magnetic field is measured. The horizontal 
component of these forces can only be fully deduced from the 
observations close to the polarity inversion line. We find that 
this component is opposed to the observed horizontal motions. 
Therefore, these Lorentz forces are a reaction of pressure- 
driven horizontal plasma motions (note the smallness of iner- 
tial effects), and correspond to magnetic free-energy storage. 
Rough estimates based on the observed parameters give insuf- 
ficient energy storage to account for the energy released in the 
flares. 

However, the rate of energy buildup was high only at the 
location of the flare activity and just preceded the start of this 
activity. Thus, the weak observed energy buildup may be an 
indication of a stronger energy buildup. It is hard to explain 


w 



Fig. 8. — Sketch showing the mass flows indicated by the observations. 
Material emerges to the surface at the location indicated with the dot in a 
circle, flows horizontally around the spot Nl, and submerges at the location 
near P3 that is indicated with a cross in a circle. The sunspot PI also moves 
horizontally but with a smaller displacement. The horizontal velocity is 
maximum close to the polarity inversion line between Nl and PI, where large 
shear of the field occurs. 


how both flares could arise from the same magnetic process, 
unless the observed vertical current’s magnitude and distribu- 
tion are not critical factors. This conclusion contrasts with 
most flare scenarios in which the vertical current plays a 
central role. 

Horizontal currents are required by MHD theory because of 
the solenoidal condition of the electric current. Horizontal cur- 
rents cannot be derived from available observations of photo- 
spheric fields, but their structure is critical for magnetic 
free-energy storage. 

A probable solution for the puzzle of the free-energy storage 
is that it occurs mostly in the form of horizontal currents. 
Important energy buildup may be associated with the horizon- 
tal electric currents due to flow of material. Lorentz forces due 
to these currents have vertical and horizontal components that 
we cannot observe, but they can do the work needed. The 
formation of horizontal currents is naturally expected when 
plasma and field emerge. Also “emergence” will force 
“horizontal flows,” and plasma will “subside” elsewhere. 
Thus, it is not possible to consider independently any of these 
three components of the overall flow pattern. 

The existence of vertical currents at the photosphere implies 
that horizontal currents must exist above the layers measured. 
Flares may be produced by a rearrangement of these currents, 
without change in the vertical currents at the photosphere. 
This process may release part of the self- and mutual-energy of 
the currents. These components can be interpreted in a circuit 
analogy, but in more realistic MHD they have to be more 
rigorously defined. Fontenla (1993) provides in his Appendix B 
such rigorous definitions valid in any full MHD case. In this 
paper Fontenla also develops the theory of current sheet 
buildup in a magnetic null configuration. Our observations 
indicate a null in the flaring region. Therefore the observed 
photospheric motions are expected to produce current sheets. 
The collapse of such sheets can produce the rapid energy 
release in the observed flares. However, for modeling of the 
MHD processes taking place in these flares a two-dimensional 
approach is not realistic and a three-dimensional approach is 
needed. In the three-dimensional case horizontal electric cur- 
rents can close consistently with vertical currents at a finite 
distance, and plasma flows can be treated consistently with 
mass conservation and gravity. 

Finally, we also note that the observed Lorentz forces pose 
problems for any attempt to construct force-free extrapo- 
lations of the observed field for assessing the magnetic free- 
energy storage. These difficulties arise because we observe 
significant non-field-aligned current that do not fit into any 
force-free extrapolation. The existence of these currents is 
logical because the observations correspond to photospheric 
layers where the plasma-beta is not very small (j8 1 for the 

observed field of 800 G at the polarity inversion line where the 
two flares occurred). 
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APPENDIX A 


CALCULATION OF THE POTENTIAL FIELD 


In the following we will assume a potential field above the photosphere. This hypothesis is known to fail in some regions because 
vector magnetograms indicate substantial vertical currents. However, measured vertical currents are often localized, and there are 
extended regions where the vertical currents are below the measurement noise level. Thus, in many cases the actual structure of the 
coronal magnetic fields may not differ too much from a potential field. 

For computing the magnetic field from the vector magnetic fields measured at the photospheric boundary we developed a scheme 
that takes into account not only the vertical, but also the horizontal components of the field. The use of the horizontal components 
partially compensates for the limited field of view and the lack of side and top constraints on the field. In theory, we would be able to 
define a unique potential function, that accurately describes the magnetic field in the current-free domain, by solving the 
equations 

= 0 , n • V<D = II • 5, , 

where indicates the measured value of the field over the surface s that fully encloses the simply connected domain. 

In practice, however, the magnetograph field of view represents only a limited portion of the enclosing surface. The observed 
rectangular portion of the photosphere may be considered as the lower boundary of a parallelepiped, in which the normal 
component of the field is unknown at the sides and top. The potential field results from solving the equations 

V^d) = 0 , e,- V<Dp = e, ‘ B^(x, y, z = 0) , 

where Bp is the vector magnetic field measured at the photospheric level, that is, the lower boundary. These equations do not fully 
determine the potential field. The additional constraints at the sides and top of the domain are necessary to discriminate among the 
set of potential fields that satisfy the previous equations. 

There have been several methods developed for computing the potential field from observations. Altschuler et al. (1977) assume 
that the normal component of the field is known at a closed spherical surface at the photosphere; a condition that the field is radial 
is imposed at an outer, “ source,” concentric spherical surface. This method resorts to a strong assumption regarding the radius of 
the “source” surface, and the fact that the field is purely radial there. Furthermore, measurements of the vertical field across the 
entire solar surface cannot be made at present, and this field is only inferred from synoptic observations over a full solar rotation. 
This method is well suited for studying the slowly varying global solar magnetic field configuration, but can hardly be justified for 
studying single dynamic active regions. Other methods for deriving the potential field resort to Fourier expansions; they assume 
periodic conditions at the sides of the domain. These conditions can hardly be justified because active regions do not repeat 
themselves. Usually the Fourier expansion methods assume that the field decreases exponentially at the upper boundary. This is 
also a very strong constraint that is not satisfied by a simple dipolar field, the simplest field representing many observed situations. 

The method that we propose consists in using the observed transverse components in selecting between all the possible potential 
fields that satisfy the previous equations. For this we construct the potential field from a distribution of arbitrarily oriented dipoles 
that are located at a constant depth below the plane of the observations. This depth, H, is chosen to be sufficiently larger than the 
separation between observed pixels, in order to minimize the nonsmoothness that arise due to the discrete representation of the 
sources of the field. The magnetic field in the domain is given by 


y, z) = X 


3r,iDrn)-D,rf 


where 


r, = {x- y - y„ z - H) 

and where the index i denotes each of the dipoles that are assumed to be the sources of the field. In our present calculations these 
dipoles are located beneath the center of each of the pixels (or rather of the group of observed pixels since we have performed some 
averaging). For our current calculations we have set H to be twice the separation between pixels, that is, H 22" and we find that 
this provides sufficiently smooth solutions to prevent sharp changes of the field in between the centers of the pixels. The three 
components of the dipole strengths are determined from the vector magnetic field data, but the horizontal components of the field 
are substantially noisier than the vertical. Moreover the observations of the horizontal field have an inherent ambiguity of 180°. 

Therefore, we performed an iterative procedure in which we initially assume that the dipoles have only vertical components. By 
inverting a linear system of equations we derive the magnitude of the dipoles such that the potential field would match the observed. 
This corresponds to solving 


Bzj = I Dz? 


- r 


r?. 


hi 


Then, the vertical component of the source dipoles that was previously determined is used to compute the horizontal components of 
the potential field corresponding to the previous iteration. This horizontal potential field is compared with the observed. At each 
iteration the 180° ambiguity in the direction of the measured field is resolved by minimizing the angle between the observed and 
computed fields. The errors of the computed horizontal field are corrected by introducing horizontal components of the source 
dipoles. These components are determined such that they, together with the previously determined vertical components, produce a 
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horizontal field that matches the observed. This corresponds to solve the following linear system of equations 


Bxj-Y.Dzr^ 

i 

Byj-Y.Dzr^ 

i 


^ = I Dx? + X Dyl , 


Since the now determined horizontal components of the dipoles also affect the vertical field, the current iteration vertical component 
of the dipoles is found by solving the following equations 


Bzj 


iDxl 


3Hx, 


''Ij 




3H^ 




'i.J 


The resulting three components of the source dipoles provide a potential field which accurately matches the observed vertical field, 
but whose horizontal field does not accurately match the horizontal fields. 

After a few iterations in this procedure the mismatch in the horizontal fields decreases substantially, then remain unchanged. 
Most of our computations use a field of view of -232" x 174" represented by 20 x 15 pixels (corresponding to 80 x 60 pixels in the 
unsmoothed data). These computations were carried in a VAX computer and require considerable memory. Calculations were also 
carried using half of the previous smoothing, and with 40 x 30 pixels. The higher resolution calculations show more detail but no 
significant differences for the selected active region. . , j u 

Our procedure was carefully tested by setting up artificial data derived for the field from a simple dipole source buried somewhere 
under the observation’s plane. When the dipole source is located too far down compared with the size of the field of view, large 
errors in the estimate of the potential field occur at some height above the plane of observations. When the source of the field is 
brought near the plane where we set up the dipoles for approximating the potential field, the artificial distribution of dipoles reduce 
to having only one dipole with the magnitude and orientation corresponding to the source of the actual field. In this case the 
approximation to the potential field converged to the real potential field with high accuracy. In all cases, when modeling the 
potential field with null error artificial data we obtained convergence to a reasonably approximate solution. We have modeled some 
artificial data in which electric currents were built in. The iterative procedure initially converges, but then it reaches a point where it 
ceases to converge. The residual rms departures of the estimated from the measured horizontal field depend on the currents; this 
residual is a measure of nonpotential behavior. In all calculations using artificial data we found that our method gave an overall 
structure of the field that behaved more like the real field than calculations using other methods, such as Fourier expansion, 
monopole distributions, or purely vertical dipole distributions. 

We conclude from these test calculations that an approximate potential can be set up by using the transverse components ol the 
measured field. This procedure selects among the potential fields that match the measured longitudinal fields over the magneto- 
graph field of view. The selection is not arbitrary but is based on the observed transverse fields. The potential field from our method 
can be obtained in all cases, has an overall structure that is reasonable, and, unlike other methods, is not much affected by the flux 
imbalance in the field of view. 


APPENDIX B 


THE EVOLUTION OF ELECTRIC CURRENTS 

We derive the electric current density from magnetograms by using Ampere’s law in its integral form, HqI =^Bdl. The 
integration we use interpolates the vector transverse fields between each group of four adjacent pixels and assumes that the current 
is uniformly distributed within the enclosed area. This gives the vertical current density values = IJS, where S is the area of the 
enclosed surface, and the current density value is assigned to the center of the area. This determination gives only a lower bound on 
the current density, because a much larger current density may be highly localized within the pixel separation and produce the same 
result. However, the net unbalanced current within the resolution area is accurately derived from our estimates (except for errors in 

the transverse field measurement). r .. j- .• t ,u 

The determination of the electric current must resort to some scheme to resolve the 180 ambiguity of the direction of the 
transverse field. Our scheme assigns the direction of the observed transverse field in such a way that the angle between observed and 
potential fields is less than 90°. This criterion is somewhat arbitrary and implies an assumption on the upper bound of the 
magnitude of the electric current. Moreover, in highly sheared regions the application of the criterion produces a definite answer 
even in the case of an ambiguity between angles of 89° or -91°. In this last case, it is hard to make any local argument for selecting 
the lower absolute value of the angle. Another scheme to resolve the ambiguity, used by Gary & Demoulin (1993) for similar data, 
resulted in basically the same electric current pattern (see their Fig. 2). This and the fact that similar patterns are observed when 
comparing different data sets gives us confidence in the significance of the patterns depicted by our Figure 5. 

Assuming an uncertainty in the transverse field of AB « 100 G, the maximum noise level in our data can be estimated to be about 
(2AB,)lano) » 10 mA m'^. However, the observations show that the actual noise level is ~4mAm 
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